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Discovery of a parsec-scale bipolar nebula around MWC349A 
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ABSTRACT 



We report the discovery of a bipolar nebula around the peculiar emission-line star MWC 349 A using archival Spitzer Space Telescope 
24 /jm data. The nebula extends over several arcminutes (up to 5 pc) and has the same orientation and geometry as the well-known 
subarcsecond-scale (~ 400 times smaller) bipolar radio nebula associated with this star. We discuss the physical relationship be- 
tween MWC 349A and the nearby BO III star MWC 349B and propose that both stars were members of a hierarchical triple system, 
which was ejected from the core of the Cyg 0B2 association several Myr ago and recently was dissolved into a binary system (now 
MWC349A) and a single unbound star (MWC349B). Our proposal implies that MWC349A is an evolved massive star (likely a 
luminous blue variable) in a binary system with a low-mass star. A possible origin of the bipolar nebula around MWC 349A is 
discussed. 
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1. Introduction 

MWC 349A is a curious star located on the sky in the direc- 
tion of the Cyg OB2 association. It was identified as a peculiar 
emission-line star in 1932 by Merrill et al. (1932), who also 
noted the presence of a second star, MWC 349B, w 2 arcsec 
west. Since that time MWC 349A attracted wide attention and 
although many interesting details have been revealed about it, 
the nature of this star remains puzzling. 

MWC 349A is one of the brightest radio stars on the sky 
(Braes et al. 1972). Its radio emission was resolved with the Very 
Large Array (VLA) in a sub-arcsecond square nebula by Cohen 
et al. ( 1 985), who suggested that the shape of the nebula might be 
due to a biconical outflow of ionized gas pinched at the waist by 
a disk seen edge-on and oriented due east-west. This suggestion 
was confirmed by White & Becker (1985), whose higher reso- 
lution and higher frequency data clearly revealed a hourglass- 
shaped bipolar nebula with a symmetry axis almost in the north- 
south direction (see also Tafoya et al. 2004). MWC 349A also 
known as a source of hydrogen (sub)millimeter recombina- 
tion lines that show maser (Martin-Pintado et al. 1989) and 
laser action (Strelnitski et al. 1996). With increasing frequency, 
these lines attain a more and more pronounced two-peaked pro- 
file that was interpreted as originating from a disk or ring- 
shaped structure in Keplerian rotation around the central star of 
mass of ~ 30 Mo (Thum et al. 1992, Ponomarev et al. 1994). 
(Sub)millimeter interferometry showed the emission to be con- 
sistent with a disk (Planesas et al. 1992), but revealed a velocity 
structure not consistent with Keplerian, which prohibits a mass 
estimate for the central object (Weintroub et al. 2008). 

Cohen et al. (1985) found a signature of interaction between 
the nebula and MWC 349B and suggested that this star forms a 
binary system with MWC 349A (see also Tafoya et al. 2004). 



Cohen et al. (1985) were also able to separate the spectra of 
the two stars and classified MWC349B as a BO III star The 
spectrum of MWC 349A is dominated by a multitude of very 
strong emission lines, which completely hide the photospheric 
lines (Cohen et al. 1985; Andrillat et al. 1996) and make a spec- 
tral classification of the star impossible. Using the spectral type 
of MWC 349B, Cohen et al. (1985) estimated a distance to the 
putative binary system of 1 .2 kpc, which is much smaller than 
the distance to Cyg OB2 accepted at this time, i.e., 2 kpc. The 
physical relationship between MWC 349A and MWC 349B was 
questioned by Meyer et al. (2002). These authors measured po- 
larization towards the two stars and found that MWC 349B is 
more highly polarized than MWC 349A. On the other hand, they 
found that several members of Cyg OB2 near to MWC 349A 
have polarization characteristics similar to those of the inter- 
stellar polarization towards this star Based on these findings, 
Meyer et al. (2002) suggested that MWC 349A and MWC 349B 
are members of the Cyg OB2 association and that they are pro- 
jected by chance near the same line of sight, with MWC 349B lo- 
cated "behind the dusty region around MWC 349 A". Meyer et al. 
(2002), however, did not explain how to reconcile the distances 
to Cyg OB2 and MWC 349B. The problem of the physical rela- 
tionship between MWC 349A and MWC 349B is very important 
for understanding the nature of the former star If both stars form 
a binary system, then MWC 349A should be a massive evolved 
star, while the distance to the system should be equal to the spec- 
troscopic distance to MWC 349B. In Sect. [3] we argue that both 
stars could indeed reside in the Cyg OB association and at the 
same time they could be physically related to each other (at least 
in the recent past). 

The bipolar radio nebula associated with MWC 349A is pro- 
duced by a constant velocity outflow of ionized gas (OInon 
1975; Hartmann et al. 1980). Widths of the [Nii] A65U emis- 
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sion line (Hartmann et al. 1980) and the hydrogen radio re- 
combination lines (Altenhoff et al. 1981) suggest that the ve- 
locity of the outflow is ~ 50 kms"'. Hartmann et al. (1980) 
argued that the low-velocity nature of the wind could be un- 
derstood if the star has a low effective surface gravity (i.e. the 
star is inflated) and interpreted MWC 349A as a supergiant star 
of PCygni type. Radio continuum observations indicate an out- 
flow mass-loss rate of ~ 10"^ Moyr"' (Hartmann et al. 1980; 
Altenhoff' et al. 1981; Dreher & Welch 1983). With this mass- 
loss rate and the low wind velocity the line emission in the opti- 
cal and the near-infrared must become optically thick if the out- 
flow originates at a radius of $ 1 AU (Altenhoff et al. 1981). 
Since the emission lines appear to be optically thin (Thompson 
et al. 1977), it was suggested that the outflow is fed by a cir- 
cumstellar disk, so that the velocity of the wind is rather set by 
the escape velocity from the disk than by that from the stel- 
lar surface (Altenhoff et al. 1981). The presence of a (nearly 
edge-on) circumstellar disk also allows to explain the coexis- 
tence of low- and high-excitation emission lines in the spectrum 
of MWC349A (Hamann & Simon 1986), the double-peaked 
profiles of optical, infrared (Hamann & Simon 1986, 1988) and 
masing (Planesas et al. 1992; Thum et al. 1992; Gordon 1992) 
emission lines, and the high level of polarization in this star 
(Elvius 1974; Yudin 1996; Meyer et al. 2002). The existence 
of the disk was confirmed with infrared speckle interferometry 
(Leinert 1986; Danchi et al. 2001; Hofmann et al. 2002), which 
revealed a flattened structure (with a major axis of ^ 60 mas or 
X 85 AU at our adopted distance of 1.4 kpc; see Sect. l3.1b lying 
in the east-west plane, i.e. oriented perpendicular to the symme- 
try axis of the bipolar radio nebula. It is widely accepted that the 
disk plays a crucial role in formation and shaping the radio neb- 
ula, whose extent along the polar axis (as seen in the VLA 2 cm 
image; see Sect.|221) is ~ 1000 AU or ^ 0.005 pc. 

In this paper, we report the discovery of a parsec-scale mid- 
infrared bipolar nebula around MWC 349 A, whose orientation 
and geometry are similar to those of the sub-arcsecond radio 
nebula. The close similarity of the morphologies of the two neb- 
ulae suggests that they are shaped by the same agent - the disk 
around MWC 349A. The newly-discovered nebula is presented 
in Sect.|2] In Sect.[3]we discuss the location and the evolutionary 
status of MWC 349 A, the possible origin of the bipolar outflow 
associated with this star, and the similarity between this outflow 
and bipolar nebulae produced by other massive evolved stars. 
We summarize in Sect.|4] 



2. Mid-infrared, radio and Ha nebulae around 
l\flWC349A 

2.1. Spitzer data 

The large-scale infrared nebula around MWC 349 was detected 
in archival data originating from observations of the Spitzer 
Space Telescope, namely within the framework of the Cygnus- 
X Spitzer Legacy Survey (Hora et al. 2008)0. This survey cov- 
ers 24 square degrees in Cygnus X, one of the most massive 
star-forming complexes in the Milky Way (e.g., Piddington & 
Minnett 1952; Reipurth & Schneider 2008), and provides im- 
ages at 24 and 70yum obtained with the Multiband Imaging 
Photometer for Spitzer (MIPS; Rieke et al. 2004) and at 3.6, 
4.5, 5.8, and 8.0yum obtained with the Infrared Array Camera 
(IRAC; Fazio et al. 2004). The resolution of the MIPS 24 and 
70jum images is a; 6 and 18 arcsec, respectively, while that of 



the IRAC images is a; 1 arcsec. Inspection of the data from this 
survey have already led to the discovery of a new Wolf-Rayet 
star (through the detection of a circular 24 fim shell and follow- 
up spectroscopy of its central star; Gvaramadze et al. 2009) and 
a second (concentric) shell around the already known infrared 
shell (Trams et al. 1998; Egan et al. 2002) surrounding the can- 
didate luminous blue variable (cLBV) star GAL 079.29+00.46 
(Gvaramadze et al. 2010a; Jimenez-Esteban et al. 2010; Kraemer 
et al. 2010). The MIPS and IRAC post-basic calibrated data 
on MWC 349A were retrieved from the NASA/IPAC Infrared 
Science Archivq3- Like many other circumstellar nebulae dis- 
covered with Spitzer (e.g. Gvaramadze et al. 2010a; Wachter 
et al. 2010; Mizuno et al. 2010), the nebula associated with 
MWC 349 A is best visible at 24-fim. It can also be seen at 70 //m, 
but in none of the IRAC images. 

The left and the middle panels of Fig.[T]present the MIPS 24 
and 70/im images of the nebula and its central star MWC 349 A 
(a white dot at the centre of the star in the MIPS 24 fim image is 
due to saturation effect). At lA^m. the nebula has a pronounced 
bipolar (hourglass-like) structure stretched almost in the north- 
south direction. The waist of the nebula is surrounded by a bright 
belt with a diameter of a; 2.7 arcmin. The nebula is somewhat 
asymmetric with respect to the polar axis. From the west side the 
lobes are outlined by almost straight filaments, with the north- 
ern and the southern ones extended, respectively, out to a 5.5 
and 6 arcmin from the waist. From the east side, the lobes are 
bounded by somewhat curved filaments. The northern filament 
extends out to a 3 arcmin from the waist and then merges with 
the bright background, while the southern one can be traced up 
to = 7 arcmin. The curvature of the latter filament suggests that 
the southern lobe might be a closed structure. 

At a 2 arcmin from the waist the eastern filament of the 
southern lobe splits into two parts and becomes concave towards 
the interior of the lobe, which might indicate an interaction be- 
tween the laterally expanding nebula and the ambient medium. 
At our adopted distance of 1 .4 kpc (see Sect. 13. lb . 1 arcmin cor- 
responds to » 0.4 pc, so that the diameter of the waist is a 1.1 
pc, while the lobes stretch out to k 2.4 - 2.8 pc from the equa- 
torial belt. Interestingly, MWC 349A is offset by a 20 arcsec (or 
K 0.14 pc) from the geometric centre of the waist, being closer 
to its eastern edge. This displacement could be understood if the 
nebula impinges on a more dense ambient medium in the east- 
ern direction. The presence of the dense material on the eastern 
side of the nebula might be inferred from the 70jum and the 17 
cm radio wavelength (see Sect. 12. 2b images of the field around 
the nebula (see the middle and the right panels of Fig.[T] respec- 
tively). The 70/Ym image also shows obvious counterparts to the 
equatorial belt and the eastern boundary of the southern lobe. 

2.2. VLA data 

To produce a long wavelength (L-band, a 17 cm) radio image 
of MWC 349 A, we used two archival datasets obtained with the 
NRAO Very Large Array (VLAJj. The datasets were acquired 
and calibrated by R. Perley in the course of a long range monitor- 
ing campaign. The data was taken on 2008 August 3 1 in the most 
compact (D) configuration at 1665.9 and 1875. 1 MHz, each with 
one intermediate frequency band with bandwidth 12.5 MHz. The 



http://www.cfa.harvard.edu/cygnusX 



" http://irsa.ipac.caltech.edu/ 

' The VLA is operated by the National Radio Astronomy 
Observatory (NRAO). The NRAO is a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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Fig. 1. Left: MIPS 24-fim image of the arcminute-scale bipolar nebula and its central star MWC 349A (a bright source in the waist 
of the nebula; note that a white dot at the centre of this source is due to saturation effect). Middle and Right: MIPS 70yum and VLA 
17 cm images of the same field. The 47" FWHM restoring beam is represented in the bottom right corner of the radio image. 



data were processed further with programs of the Astronomical 
Image Processing System (AIPS) in the usual manner. First, 
the two datesets were combined (using DBCON) and then im- 
aged with IMAGR. One cycle of self caUbration (with CALIB) 
improved the image quality. The image, part of which is pre- 
sented in Fig.[I] was restored with a circular Gaussian beam with 
FWHM 47". There are no clear indications of the presence of ra- 
dio emission associated with the mid-infrared nebula. 

We also produced a radio image of MWC 349A with higher 
resolution taken in the most extended (A) configuration at a 
shorter wavelength (U-band, 2 cm) from archival data taken 
on 1983 October 30 at a frequency of 14.94 GHz. The obser- 
vations are reported and discussed by White & Becker (1985) 
and we refer to this publication for details on the data and 
their processing. Additional use of self-calibration resulted in 
an improved image. This image (presented in Fig.|2) shows 
the well-known subarcsecond-scale radio nebula associated with 
MWC 349A. Comparison of Fig.|2]with the MIPS 24 //m image 
of the arcminute-scale nebula shows that both nebulae have the 
same orientation and bipolar morphology. Note that by coinci- 
dence the angular size of the field shown in Fig. |2] turns out to 
be equal to that of the saturated pixel (the white dot) in the very 
centre of the 24 /im image of MWC 349 A. 



2.3. Ha data 

The parsec-scale nebula around MWC 349A was also detected 
in the Ho- narrow -band CCD imaging survey of B[e] stars by 
Marston & McCollum (2008). These authors noted that "an ap- 
parent thin shell of approximately 2 '5 diameter exists around 
MWC 349A to the noith". Comparison of their Fig. 4a with the 
24yL(m image in Fig. [T] shows that a thin Ha filament west of 
MWC 349A coincides with the western edge of the equatorial 
belt, while an arc-like diffuse Ha emission to the east shows a 
good correspondence with the eastern edge of the waist. 
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Fig. 2. VLA 2 cm image of the subarcsecond bipolar nebula 
around MWC 349A (indicated by a cross). The angular size of 
the field shown in this figure is equal to that of a white spot in 
the centre of the 24-fim image of MWC 349 A (see the left panel 
of Fig.[T]). The 0" 1 1 FWHM restoring beam is represented in the 
bottom left corner. 



3. Discussion 

Above we have reported the discovery of a parsec-scale 
mid-infrared nebula around the enigmatic emission-line star 
MWC 349 A, whose hourglass-like structure and orientation are 
similar to those of the well-known subarcsecond radio nebula 
associated with this star. The morphology of the nebula is very 
similar to that of some planetary nebulae (e.g. the Engraved 
Hourglass Nebula; Sahai et al. 1999), nebulae around symbi- 
otic systems (e.g. the Southern Crab Nebula; Corradi & Schwarz 
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1993), and nebulae associated with young stellar objects (e.g. 
the disk/outflow radio source I in the Kleinmann-Low Nebula in 
Orion; Plambeck et al. 2009; Matthews et al. 2010) and some 
evolved massive stars [e.g. the cLBV Sher25 (Brandner et al 
1997) and the famous LBV //Carinae]. Correspondingly, three 
main interpretations of the subarcsecond bipolar nebula asso- 
ciated with MWC349A were proposed: (i) planetary nebula 
(Ciatti et al. 1974), (ii) photo-evaporation induced outflow from 
an accretion disk around a massive pre-main sequence star (e.g. 
Thompson et al. 1977), and (iii) outflow from an (excretion) disk 
around an evolved massive star (Hartmann et al. 1980). 

Hofmann et al. (2002) thoroughly analyzed these three pos- 
sibilities and came to the conclusion that MWC 349A is most 
likely an evolved massive star (a B [e] supergiant) residing in the 
Cyg OB2 association. Similar proposals on the evolutionary sta- 
tus of MWC 349A were also put forward by Swings & Struve 
(1942), Baldwin et al. (1973), Hartmann et al. (1980), Herzog et 
al. (1980), Andrillat et al. (1996), and Lamers et al. (1998). The 
membership of the Cyg OB2 association seems, however, to be 
at variance with the possibility that MWC 349A forms a physical 
system with MWC 349B (Cohen et al. 1985; Tafoya et al. 2004), 
since the spectroscopic distance to the latter star of ~ 1 .2 kpc is 
significantly smaller than the today's widely accepted distance 
to Cyg OB2 of 1 .7 kpc (Knodlseder 2000; Reipurth & Schneider 
2008). Before discussing the possible origin of the bipolar neb- 
ula associated with MWC 349 A, we have to understand whether 
the physical relationship between the two stars could be consis- 
tent with their membership of the Cyg OB2 association. 

3.1. Birth cluster and the evolutionary status of MWC 349 A 

Let us assume that MWC 349A is a massive evolved star as ar- 
gued by Hofmann et al. (2002). In this case, the question arises: 
where is the birth cluster of this star? Indeed, there is grow- 
ing acceptance that most (or probably all) massive stars are 
formed in the clustered way (Lada & Lada 2003) and subse- 
quently found themselves in the field either because of dynami- 
cal few-body encounters in the parent clusters, binary-supernova 
explosions and/or due to rapid cluster dissolution (e.g., Kroupa 
& Boily 2002; de Wit et al. 2005; Schilbach & Roser 2008; 
Gvaramadze & Bomans 2008b; Pflamm-Altenburg & Kroupa 
2010; Gvaramadze et al. 2010b; Weidner et al. 201 1 fl One can 
therefore expect that the parent cluster of MWC 349A should 
be nearby, unless this star is a high-velocity runaway. Since 
MWC 349A is located in the direction of one of the most com- 
pact and massive associations in the Milky Way - the Cyg OB2 
association (Knodlseder 2000), it is natural to assume that it is a 
member of this stellar system. 

If MWC 349A is indeed a member of Cyg OB2 then the 
small angular separation between MWC 349A and MWC 349B 
might simply be due to a chance projection. The line of sight 
towards MWC 349A is nearly tangential to the local (Orion) 
spiral arm (whose extent in this direction is ~ 4 - 6 kpc; e.g. 
Russeil 2003), which makes the chance projection quite pos- 
sible. The SIMBAD databas^ lists three O stars within 10 ar- 
cmin from MWC 349 A, one of which, the 08 III (Negueruela et 
al. 2008) star 2MASS J20323843+4040445, is sepai-ated from 
MWC 349A by only « 1 .8 arcmin and is projected on the west- 
ern edge of the 24 //m nebula (not far from the equatorial belt). 
We note, however, that recent very accurate parallax measure- 



'' See Parker & Goodwin (2007), Seller, Heydari-Malayeri & 
Gouliermis (201 1), and references therein for a different point of view. 
^ http://simbad.u-strasbg.fr/simbad/ 



ments of five massive star-forming regions towards CygnusX 
showed that four of them are located at a distance of 1 .4 + 0. 1 kpc 
(Rygl et al. 2012). At the same time, there are strong indications 
that most molecular clouds in Cygnus X form a coherent com- 
plex located at the same distance as Cyg OB2 (Schneider et al. 
2006), which implies that Cyg OB2 is at the distance of 1 .4 + 0. 1 
kpc as well (cf. Hanson 2003). From this it follows that the phys- 
ical relationship between MWC 349A and MWC 349B and their 
membership of Cyg OB2 do not contradict each other, which in 
turn implies that MWC 349A should be as old as MWC 349B 
(i.e. ~ 5 Myr). The obvious consequence of the latter implica- 
tion is that MWC 349A cannot be a young stellar object, nor a 
post-AGB stai". 

To check whether the spectral classification of MWC 349B 
is consistent with the possibility that this star is a member of 
Cyg OB2, we use its visual magnitude and extinction of, respec- 
tively, V = 14.3 mag and Ay = 8.6 - 8.9 mag (Cohen et al. 1985) 
to derive the absolute visual magnitude My - -(5.0 -^ 5.3) mag, 
which agrees well with that expected for a BO III star. 

Similarly, assuming the visual magnitude of MWC 349A of 
14 mag (Cohen et al. 1985) and adopting Av/ - 10.0 - 10.6 ma^ 
(Cohen etal. 1985; Kelly etal. 1994), one has Mv/ = -(6.7^7.3) 
mag. The absence of He ii lines and the presence of strong He i 
emission lines in the spectrum of MWC 349A imply that the 
effective temperature of this star should be in a range from 
~ 20 000 to 28 000 K, which corresponds to a bolometric cor- 
rection of -2.5 ± 0.4 mag (Hofmann et al. 2002). Thus, one de- 
rives a the stellar luminosity of log(L/L0) - (5.6 -H 5.9) + 0.2 (cf. 
Hartmann et al. 1980; Hofmann et al. 2002). We caution that this 
luminosity estimate should be considered approximate because 
of the large {k 1 -2 mag) photometric variability of MWC 349A: 
in the B-band its magnitude varies from ^ 14 to 16 mag (Gottlieb 
& Liller 1978), while in the /J-band from ^ 9.5 to 10.5 mag 
(Jorgenson et al. 2000). On the other hand, the good correlation 
between the brightness in the V and R bands in MWC 349A (see 
Fig. 1 in Yudin 1996) implies that V = 14 mag corresponds to 
the minimum brightness of this star, so that its actual luminosity 
could be even higher 

The high luminosity of MWC 349 A, along with its signif- 
icant spectroscopic (e.g. Andrillat et al. 1996) and photomet- 
ric variability, and the presence of the parsec-scale (bipolar) 
circumstellar nebula are typical of evolved massive stars be- 
longing to the class of LBVs (e.g. Humphreys & Davidson 
1994; Bohannan 1997). We therefore concur with Hofmann et 
al. (2002) that MWC 349A is an evolved massive star, likely an 
LBV (cf. Hartmann et al. 1980). This inference is supported by 
the possible existence of a high-velocity (^ 260 km s"' ) compo- 
nent in the bipolar outflow (Tanaka et al. 1985), which might cor- 
respond to the LBV wind concentrated along the polar axis (we 
will return to this point in Sect. l3.3b . Assuming that MWC 349A 
is an LBV of age of 5 Myr and using the above estimate of its 
luminosity, one can infer the zero-age main-sequence mass of 
this star of ~ 40 Mq. 

3.2. MWC349A as a runaway 

MWC349A is separated from the geometric centre of the 
CygOB2 association by ^ 24 pc in projection. CygOB2 con- 
tains ~ 100 O stars or stars with O-type progenitors and its half 
light radius is ~ 6 pc (Knodlseder 2000). The age of the associ- 
ation is ~ 5 Myr (see Gvaramadze & Bomans 2008a and refer- 



* Note that the interstellar extinction towards the Cyg 0B2 associa- 
tion is very patchy and ranges from ~ 5 to 20 mag (Knodlseder 2000). 
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ences therein). Assuming that the association expands with a ve- 
locity equal to its velocity dispersion (^ 2.4 kms"'; Kiminki et 
al. 2007), one finds that the majority of massive stars in Cyg OB2 
were originally concentrated in a region of radius of < 1 pc. It 
is therefore plausible that several Myr ago the stellar number 
density in the core of Cyg OB2 was high enough to ensure that 
close dynamical encounters between its members were frequent, 
which is the necessary condition for effective production of run- 
aways. Currently only one runaway associated with CygOB2 
is known - the 04If star BD+43° 3654 (Comeron & PasquaH 
2007; Gvaramadze & Bomans 2008a). Let us check whether 
MWC 349A might be another example of a massive star running 
away from Cyg OB2. 

To check this possibility, we use the proper motion measure- 
ments for MWC 349 A given in Rodriguez et al. (2007). To con- 
vert the observed proper motion (fia cos 6 - -3.1 + 0.5 mas 
yr"' and fig - -5.3 + 0.5 mas yr"') into the transverse pecu- 
liar velocity, we use the Galactic constants /?() = 8.0 kpc and 
00 = 240 km s"' (Reid et al. 2009) and the solar peculiar mo- 
tion (Ue,Ve,We) = (11.1, 12.2,7.3)kms-' (Schonrich et al. 
2010). The derived velocity components in Galactic coordinates 
are V; = -7.3+6.8 kms"' andv/, = 3.0+3.3 kms"' (for the error 
calculation, only the errors of the proper motion measurements 
and the distance estimate were considered). With the transverse 
peculiar velocity of iv ~ 8 km s"', MWC 349 A would need ^ 3 
Myr to travel from the centre of the association to its current 
position. 

Figure|3] shows that the peculiar velocity of MWC 349 A has 
somewhat "wrong" orientation, i.e. it does not point away from 
the centre of the CygOB2 association. To explain this mis- 
alignment, one can assume that Cyg OB2 has a peculiar (trans- 
verse) velocity of ~ 5-10 kms"' in the northwest direction 
(cf. Hoogerwerf et al. 2001). The peculiar velocity of just this 
magnitude and orientation is required to explain the relative po- 
sition of the runaway star BD+43° 3654 and two young pul- 
sars on the sky (Gvaramadze & Bomans 2008a). Peculiar ve- 
locities of ~ 10 kms ' are typical of the OB associations near 
the Sun (de Zeeuw et al. 1999). More importantly, three star- 
forming regions adjacent to the CygOB2 association (DR20, 
DR21 and IRAS 20290+4052) have peculiar (transverse) veloc- 
ities of * 6-9 kms ' directed northwest and west (Rygl et al. 
2012; see also Fig.[3]l. 

To the peculiar transverse velocity of MWC 349 A one should 
add the peculiar radial velocity, Vjad, which can be derived from 
the recession velocity of this star of 8 - 12 kms"' (e.g. Thum 
et al. 1992; Gordon et al. 2001; Gordon 2003). After correc- 
tion for differential Galactic rotation and solar peculiar motion, 
one has v'rad ~ 18-22 kms"' and the total space velocity 
Vtot ~ 20 - 23 kms"'. Although this velocity does not meet the 
usual criterion to define runaway stars (i.e. Vtot > 30 kms"'; 
Blaauw 1961), we note that the ejection of stars with low veloc- 
ities is common (e.g. Gies 1987; Kroupa 1998; Gvaramadze & 
Bomans 2008b), so that any star unbound from the parent cluster 
should be considered as a runaway, independently of its peculiar 
velocity. 

On the other hand, the derived total velocity of MWC 349A 
is too large to be consistent with the observed linear separa- 
tion between MWC 349A and MWC 349B {^ 3300 AU) and the 
possibility that these stars form a bound system. This inference 
follows from the results of numerical scattering experiments by 
Kroupa (1998), who showed that a wide binary system can sur- 
vive the ejection process only if its space velocity is less than 
the orbital velocity. Simple estimates show that for any reason- 
able masses of MWC 349A and MWC 349B the orbital velocity 
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Fig. 3. Midcourse Space Experiment (MSX) 8.3 yum image of 
the Cyg OB2 association and its environments. The approximate 
boundary of the association is shown by a dashed circle of a 
diameter of a; 2°. The position of MWC 349 A is marked by a 
diamond, while that of the nearby cLBV GAL 079.29+00.46 
is indicated by a circle. The arrows show the space motions 
of MWC 349 A and three star- forming regions (indicated by 
squares). See text for details. 



of the binary is an order of magnitude smaller than its space ve- 
locity. A single exception (i.e. a high-velocity wide binary) that 
occurred in the experiments by Kroupa (1998) was interpreted as 
the outcome of "a complex high-order interaction, that resulted 
in two stars being ejected on essentially identical trajectories". 
Another possibility to produce a wide runaway binary is through 
the ejection of a compact, initially stable hierarchical triple sys- 
tem, which in the course of evolution of its components becomes 
unstable and dissolves into a binary star and a single unbound 
star We further discuss this possibility in Sect. 13. 3 1 because it 
might be directly related to the formation of the bipolar nebula 
around MWC 349A. 

To conclude, we note that MWC 349A is located at only 
R! 20 arcmin (or ^ 8 pc in projection) from the cLBV 
GAL 079.29+00.46 (see Fig.©. The distance estimates for 
GAL 079.29+00.46 (Voors et al. 2000) suggest that this stai- 
could be a member of CygOB2 as well. GAL 079.29+00.46 
is surrounded by a curious circumstellar nebula consisting of 
two concentric circular shells (see, e.g.. Fig. 2j in Gvaramadze 
et al. 2010aj3- The diameter of the inner (bright) shell is x; 1.4 
pc (at the assumed distance to GAL 079.29+00.46 of 1.4 kpc), 
while that of the outer shell is x 2.8 pc. We speculate that the 
nebula around GAL 079.29+00.46 might actually have a bipo- 
lar geometry (similar to that around MWC 349 A) with the po- 
lar axis parallel to our line of sight. In this case, the inner shell 
would correspond to the circular waist of the nebula around 
MWC 349 A (whose diameter is ^ 1.1 pc), while the outer one to 



' For another example of a two-shell circumstellar nebula produced 
by the cLBV MN112 see Gvaramadze et al. (2010c). 
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Fig. 4. MIPS 24yum (left) and IRAC 8yL/m images of the bipo- 
lar nebula around the candidate LBV MN13 (Gvaramadze et al. 
2010a). The orientation and the scale of the images are the same. 



the bipolar lobes (whose diameter is x; 2.4 pc). If subsequent 
studies of MWC 349A will prove its LBV nature, then along 
with GAL 079.29+00.46 it would represent one more example 
of massive stellar twinning (Walborn & Fitzpatrick 2000). Two 
other well-known examples of such twinning are the pairs of 
(c)LBVs: HD 168607 and HD 168625 (separated by only ^ 1.6 
arcmin) and AG Car and Hen 3-519 (separated by =s 16.2 ar- 
cmin). 

3.3. Possible origin of the bipolar nebuia around MWC349A 



In Sect.[37T]we suggested that MWC 349A is likely an LBV star 
Observations show that almost all LB Vs are surrounded by com- 
pact nebulae with linear extent ranging from ~ 0.1 to several 
pc (e.g. Weis 2001; Clark et al. 2005). These nebulae display 
a wide diversity of shapes, from circular to bipolar and triple- 
ring forms (e.g. Nota et al. 1995; Smith 2007; Gvaramadze et 
al. 2010a). As a rule, the youngest LBV nebulae have bipo- 
lar morphologjo (see Table 1 in Weis 201 1). Three well-known 
examples of young (^ 10'* yr) bipolar nebulae associated with 
Galactic (c)LBVs are the Homunculus nebula of the LBV ?/Car 
(e.g. Morse et al. 1998), the two-lobe nebula around the cLBV 
Sher25 (Brandneret al. 1997), and the triple-ring system around 
the cLBV HD 168625 (Smith 2007). Two other striking exam- 
ples of hourglass-like nebulae produced by evolved massive 
stars were recently discovered with Spitzer (Gvaramadze et al. 
2010a). One of them (Fig.© is created by the cLBV MN13 
(Gvaramadze et al. 2010a; Wachter et al. 2011) and the second 
one (Fig.|5]l by the blue supergiant star MN18 (Gvaramadze et 
al. 2010a). The MIPS 24 //m image of the nebula around MN13 
shows that the brightest part of this nebula is reminiscent of 
the triple -ring system produced by the cLBV HD 168625 (see 
Fig. Id in Smith 2007), while the IRAC 8/^m image shows that 
the waist of the nebula is surrounded by a belt similar to that 
of the 24 //m nebula around MWC 349 A (the possible origin of 
such belts is discussed in Sect. l3.4b . The MIPS and IRAC images 
of MN18 show that the central star is surrounded by a toroidal 
structure, which most likely is responsible for the pinching the 
waist of the nebula. 

The widely accepted explanation of the origin of two-lobe 
nebulae around evolved massive stars is that the (spherically 
symmetric) stellar wind is collimated by an aspherical (disk- 
like) circumstellar environmenlQ (e.g. Frank et al. 1995; Langer 



^ Note that such nebulae may become more spherical with time be- 
cause of the lateral expansion of their polar lobes (cf. Sect. l3.4t . 

' Another possibility is that the stellar wind is intrinsically bipolar 
(Owocki & Gayley 1997; Maeder & Desjacques 2001). 



Fig. 5. MIPS 24yum (left) and IRAC 8yum images of the bipo- 
lar nebula around the blue supergiant MN18 (Gvaramadze et al. 
2010a). The inset shows the waist of the nebula in a different in- 
tensity scale. The orientation and the scale of the images are the 
same. 



et al. 1999). Such structures could form around single massive 
stars if they reach the Q limit (i.e. if the centrifugal force and 
the radiative acceleration cancel out gravity at the stellar equa- 
tor), so that the stellar wind becomes strongly confined to the 
equatorial plane (Langer 1997, 1998). Alternatively, the origin 
of disk-like density enhancements could be caused by various bi- 
nary interaction processes, such as (i) common envelope ejection 
(e.g. Morris 1981), (ii) mass loss through the second Lagrange 
point L2 (e.g. Livio, Salzman & Shaviv 1979), and (iii) focusing 
of the primary's wind towards the orbital plane by the gravita- 
tional field of the companion star (e.g. Fabian & Hansen 1979; 
Mastrodemos & Morris 1999). All of these processes can pro- 
duce outflowing (excretion) disks in the orbital plane. 

In Sect. l3.2l we mentioned that the large linear separation be- 
tween MWC 349A and MWC 349B could be reconciled with 
their physical relationship if both stars were originally members 
of a hierarchical triple system. In support of this possibility we 
note that Jorgenson et al. (2000) reported detection of a 9-year 
periodicity in the red light variations of MWC 349 A, which im- 
plies that the star itself might be a binary system (Jorgenson et al. 
2000; Hofmann et al. 2002). Below we show how an initially sta- 
ble runaway triple system could be dissolved into a binary sys- 
tem and a single unbound star, and suggest that the gravitational 
focusing of the slow wind emanating from MWC 349A during 
the red supergiant phase might be responsible for the origin of 
disk-like circumstellar environment, which in turn is responsible 
for the origin and shaping of the observed bipolar nebula. 

Imagine that the inner binary of the triple system was com- 
posed of MWC 349A (with the zero-age main-sequence mass, 
nil, of 4OM0) and a low-mass star (with the mass, 1112, of, say, 
5 Mq), while MWC 349B (with the mass nii - 20 Mq) was in 
a wide orbit. For the sake of simplicity, we assume that ini- 
tially both orbits were circular This triple system would be hard 
enough to survive dynamical ejection from CygOB2 if the or- 
bital separation of the outer binary flout ^100 AU, while the sys- 
tem itself would be stable if the ratio of the outer to inner orbital 
periods X = Pout/An ^ 4 (Kiseleva et al. 1994), i.e. if the orbital 
separation of the inner binary flin < 35 AU. A star with initial 
mass of 40 Mq evolves through the red supergiant and yellow 
hypergiant phases (e.g. Humphreys 1991; Oudmaijeret al. 2009) 
and at the beginning of the LBV phase has a mass of ^ 20 Mq. 
Adopting this mass as the current mass of MWC 349 A, one has 
the current mass of the inner binary of 25 Mq, so that the 9-year 
orbital period of this binary corresponds to an orbital separation 
of ^ 13 AU (i.e. the secondary star orbit within the disk; cf. 
Jorgenson et al. 2000). The large orbital separation implies that 
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the inner binary did not experience the common envelope phase, 
i.e. the mass decrease of the primary star (MWC 349 A) caused 
by (spherically symmetric) stellar wind mass loss led only to the 
increase of the orbital period and separation of the binary by fac- 
tors of (45/25)" and 45/25, respectively. Thus, the initial orbital 
period and separation of the inner binary were w 3 yr and ^ 8 
AU, respectively. 

Let us assume that initially the triple was stable (i.e. X > 
4). In the process of ejection from CygOB2 the outer binary 
acquired an eccentricity gout (cf. Hoffer 1983; Hills 1975) and 
the perturbed triple system remained stable provided that 



< 1 - 



Fix) 



where 



Fix) 



2Ax^ + l.lx + 3.7 
(1 +x3)i/3(l+jc) 



(1) 



(2) 



and X = [(mi + m2)/mo,y^^ (Eggleton & Kiseleva 1995; Livio & 
Pringle 1998). Assuming X = 10 (i.e. Pout = 30 yr and flout ~ 40 
AU), one finds from Eqs. ([TJ and ^ that the system is stable 
if Cout ^ 0.42. We assume that gout obeys this condition (say, 
Cout = 0.30), so that the triple remains intact until the primary of 
the inner binary, MWC 349 A, evolved off the main sequence and 
lost a half of its initial mass during the red supergiant and yellow 
hypergiant phases. In response to the mass loss the orbital peri- 
ods of the outer and the inner binaries increased, respectively, by 
factors of (65/45)- = 2.09 and (45/25)- = 3.24 (we assume that 
1112 and 1112 are constant), so that X decreased by a factor of 0.64 
and the condition for stabiUty of the triple becomes gout ~ 0.19. 
The spherically symmetric mass loss does not affect the eccen- 
tricity of the (inner and outer) binaries (e.g. Eggleton 2005). This 
and our assumption of gout = 0.30 implies that reduction of the 
mass of MWC 349A made the system unstable. According to 
the numerical experiments by Kiseleva et al. (1994), an unstable 
triple dissolves within ~ 100 crossing times (which in our case 
corresponds to ~ 10"* yr) and leaves behind a binary system and 
a single unbound star Thus, we propose that this binary system 
corresponds to MWC 349A and its putative low-mass compan- 
ion, while MWC 349B is the third component of the triple, which 
became unbound in the recent past. 

We speculate that the slow (red supergiant) wind from 
MWC 349A was focused by the gravitational field of the 
companion star to produce a circumbinary disk (see, e.g., 
Mastrodemos & Morris 1999), which allows the current fast 
(LBV) wind of MWC 349A to expand only in the polar direc- 
tions and is responsible for the B[e] morphology of the stel- 
lar spectrum. We speculate also that the hourglass-shaped low- 
velocity outflow observed as the subarcsecond radio nebula orig- 
inates because of gas entrainment in the boundary layer between 
the stellar wind and the disk, while the wind itself is confined 
near the polar axis of the slow outflow and is mostly unobserv- 
able. The possible existence of a high-velocity (x; 260 kms"') 
component in the bipolar radio nebula (Tanaka et al. 1985) con- 
forms with our scenario and could be attributed to the stellar 
wind channelled along the polar axis. The energy injected in the 
ambient medium by the stellar wind leads to the origin of parsec- 
scale laterally expanding lobes, which we observe as the mid- 
infrared hourglass nebula. For a wind velocity of 260 km s ' and 
a linear extent of the lobes of x; 3 pc, one derives a kinematic age 
of the nebula of w lO"* yr. Note that this time-scale is comparable 
to the dissolution time of the triple system (see above), so that 
it is likely that MWC 349B became unbound from MWC 349A 
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Fig. 6. Schematic of the proposed origin of the belt around the 
waist of the 24-fim bipolar nebula associated with MWC 349 A 
(not to scale): (1) the incident hemispherical shock wave (one of 
the two lobes of the nebula); (2) the ring-like Mach shock wave 
(the belt); (3) the tangential discontinuity; (4) the reflected shock 
wave. See text for details. 



only recently. This explains why these stars are still close to each 
other. 



3.4. Origin of tine equatoriai beit of fhe parsec-scaie nebuia 

Now we discuss the possible origin of the belt around the 
waist of the mid-infrared nebula associated with MWC 349 A (cf. 
Gvaramadze 1998). Near the waist each lobe of the bipolar neb- 
ula can be considered as a hemispherical shock wave, which col- 
lide with each other in equatorial plane (see Fig.|6]for schematic 
of this interaction). During the collision, the equatorial plane acts 
as a wall from which the incident waves are reflected. Initially, 
the reflection occurs in a regular fashion: the incident and re- 
flected shock waves touch the "wall". However, when the angle 
between the wall and the tangent to the incident shock becomes 
greater than some critical value icp^- ~ 35° for a medium with an 
adiabatic index of 5/3), the regular reflection becomes impossi- 
ble (e.g. Courant & Friedrichs 1948). The incident shock wave is 
reflected without reaching the equatorial plane. In addition to the 
incident and reflected shock waves, another shock wave, which 
propagates parallel to the wall, appears (Fig.|6). Such a reflection 
is called a Mach reflection, and the emerging new shock wave (in 
our case, in the shape of a belt; cf. Chemin et al. 1995) is called 
the Mach shock wave. The three shock waves intersect along a 
single (circular) line, from which the surface of a tangential dis- 
continuity separates the medium behind the Mach shock wave 
from the medium that passes through the incident and reflected 
ones. 

The velocity of the Mach shock wave vm is related to the 
velocity of the incident shock wave vi by vm = Vi / sin 0, where </> 
is the angle between the tangent to the front of the incident shock 
and the equatorial plane of the nebula. The Mach wave reaches 
its maximum velocity at the beginning of the Mach reflection, 
i.e., at (/> = ipcr'. v™" ^ 1.7vi. As the angle (p increases further, 
I'M decreases, and the Mach wave degenerates. Using the 24 yum 
image of the nebula, we estimated that s; 50°, so that we expect 
that the belt expands with a velocity of w 1.3 times exceeding the 
local lateral expansion velocity of the bipolar lobes. 

Our interpretation of the belt around the waist as a Mach 
shock is testable. In Sect. 12. 31 we mentioned that the waist has 
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an optical counterpart and that its western edge coincides with a 
prominent Ha filament (Marston & McCollum 2008). If the belt 
is indeed the Mach shock then we expect that high-resolution 
spectroscopy of the Ha nebula would show that the expansion 
velocity of the belt is higher than that of the adjacent parts of the 
lobes. 



4. Summary 

We have reported the discovery of a parsec-scale bipolar neb- 
ula around the curious emission-line star MWC 349A. The mor- 
phology of the nebula and its orientation are similar to those of 
the well-known (~ 400 times smaller) radio nebula associated 
with MWC 349A. The close similarity of the two nebulae sug- 
gests that they are shaped by the same agent - the disk around 
MWC 349A. Significant photometric and spectroscopic variabil- 
ity of MWC 349A along with the presence of the parsec-scale 
bipolar nebula strongly argue that this star is an LBV. We have 
discussed the physical relationship between MWC 349A and the 
nearby BO III star MWC 349B and suggested that both stars were 
members of a hierarchical triple system, which was dynamically 
ejected from the core of the CygOB2 association several Myr 
ago. The stellar wind mass loss from the most massive (-40 Mq) 
star in the triple (MWC 349 A) made the system unstable, so that 
it dissolved into a binary (MWC 349 A) and a single unbound star 
(MWC 349B). The binary nature of MWC 349A is supported by 
the existence of a 9-year periodicity in the red light variations of 
this star We proposed that the gravitational focusing of the slow 
(red supergiant) wind from MWC 349A by the companion star 
led to the origin of a circumbinary disk, which currently colli- 
mates the fast (LBV) wind of MWC 349A and is responsible for 
the B [e] morphology of the stellar spectrum. 
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